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SUMMARY 


The dimensionless; transformed; nonlinear differential equation 
developed in NASA TR R -11 for describing the approximate motion and heat- 
ing during entry into planetary atmospheres for constant aerodynamic 
coefficients and vehicle shape has been modified to include entries during 
vhich the aerodynamic coefficients and the vehicle shape are varied. The 
generality of the application of the original equation to vehicles of 
arbitrary veight; size; and shape and to arbitrary atmospheres is retained. 
A closed-form solution for the motion; heating; and the variation of drag 
loading parameter m/Cj^A has been obtained for the case of constant; 
maximum resultant deceleration during nonlifting entries. This solution 
requires certain sin5>lifying assumptions vhich do not compromise the 
accuracy of the results. 

The closed-form solution has been used to determine the variation of 
m/CpA required to reduce peak decelerations and to broaden the corridor 
for nonlifting entry into the earth’s atmosphere at escape velocity. The 
attendant heating penalty is also studied. 


INTRODUCTION 


A number of studies have been made of the motion and heating of 
vehicles entering planetary atmospheres. In one such study (ref. l) an 
approximate analytical method for studying entry into any planetary 
atmosphere vas developed and applied to vehicles of arbitrary constant 
weight; size; and shape with constant aerodynamic coefficients. In 
reference 2 the method of reference 1 was applied to the study of corri- 
dor depth and guidance requirements. Also presented in reference 2 is a 
limited discussion of the effect on corridor depth of varying (modulating) 
the aerodynamic coefficients and the vehicle shape during entry. Other 
investigators have made studies of the trajectories for modulated entries. 
For example; entries at circular and supercircular velocities with varying 
lift-drag ratioS; but with constant drag; have been studied in reference 3 ; 
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entries at circular and supercircular velocities with a particular 
variation of drag vith lift have been studied in references hy 5^ and 6; 
and nonlifting entries at circular velocity (and one case of entry at 
supercircular velocity) vith varying drag have been studied in reference 7 . 
The analyses of references 3 through 7 s-ll employ time as the independent 
variable and separate sets of calculations are reqiiired for each specific 
vehicle and planetary atmosphere. 

The puip)Ose of the present investigation is to extend the method of 
reference 1 to include vehicles vith variable aerodynamic coefficients 
and variable vehicle shapes. The general applicability of the method to 
vehicles of arbitrary veight^ size^ and shape and to arbitrary atmospheres 
is retained. Calculations are made to illustrate the effect of varying 
■yshicle drag coefficient and area on the corridor depth and heating 
characteristics of nonlifting vehicles entering the earth^ s atmosphere at 
escape velocity. The type of variation considered is that vhich reduces 
the maximum resultant deceleration to some specified constant value during 
the modulation period. The variation of vehicle drag coefficient and area 
req^uired to maintain a constant resultant deceleration is determined. 
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NOTATION 

resultant deceleration^ ft sec”^ 
reference area for drag^ sq ft 

constant defined by the deceleration at the beginning of 
modulation^ 

drag coefficient ^ 2D/ pV^A 
drag force ^ lb 

conic perigee parameter (eq. (22))^ dimensionless 

gravitational acceleration^ ft sec*“^ 

gravitational conversion constant^ 32*2 ft sec"^ 

deceleration in g*s^ a/ g 

lift force j lb 

mass of vehicle ^ slugs 

convective heating rate per unit area at the stagnation 
point/ Btu ft"2 sec""i 


4T- ^ > 


dimensionless function proportional to convective heating 
rate (eq. (26)) 

total convective heat ahsorhed per unit area at the stagnation 
pointy Btu ft’^ 

dimensionless fimction proportional to total heat ahsorhed^ 
(eq. (28)) 

distance from planet center ^ ft 

radius of curvature of vehicle surface^ ft 

dimensionless radius used in the definition of 
(See eq. (23) .) 

time; sec 

tangential velocity component normal to a radius vector^ 
ft sec"^ 

dimensionless velocity ratio ; u//^ 

radial velocity component; ft sec"^ 
resultant velocity; \/u^ + = u/cos 7 

dimensionless velocity ratiO; V/^/^ 
weight of vehicle at earth* s surface; mgc^ lb 
altitude; ft 

dimensionless function of u (eq. (3)) 
atmosphere density decay parameter; ft"^ 

flight-path angle relative to the local horizontal; negative 
for descent 

vehicle parameter; m/CpyA 

corridor depth between conic perigee altitudes; statute miles 
atmosphere density; slugs ft“^ 

mean value for exponential approximation to atmosphere 
density-altitude relation; slugs ft"^ 
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initial value 
maximum value 
modulated- entry value 
surface of a planet 
overshoot houndary 
conic perigee point 
undershoot boundary 
unmodulated- entry value 

point in the trajectory where modulation begins 
point in the trajectory where modulation ends 
phase of a trajectory 


Superscripts 


differentiation with respect to u 


MALYSIS 


In reference 1 the two-component equations of motion for entry into 
planetary atmospheres 

- ^ = g - — - — 7+5 sin 

dt r m m 

+ HY = - 5 cos 7 + “ tan 7 '^ 

dt r m D J 

are transformed into a single, ordinary, nonlinear, differential equation 
with the aid of certain sinrplifying assumptions and by the introduction 
of the dimensionless independent variable 
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and the dimensionless dependent variable defined by 



(3) 


The resulting differential equation for the Z function is 

uZ* ’ - ( Z,* - — ) - ^ ~ ^ cos'^y + — cos^7 = 0 (4) 

\ u / uZ ^ D 

vhere cos 7 = \l 1 ^ sln^j can he expressed in terms of Z and Z^ through 
the equation for the flight -path angle 


sin 7 = Z^ - -i (5) 

u 

In the analysis in reference 1 the terns m/CpA and L/D are constant 
throughout the entry trajectory. 

When the teins m/CpA and L/D are allowed to vary during entry (are 
modulated) ^ the differential equation for the Z function is^ as derived 
in appendix given hy 



1 - u^ 
uZ 


cos^7 + “ cos^7 = 0 


(6) 


where 


A = m/CpA 


and 


^ d(m/Ci>A) 
du du 


(7) 


The attendant expression for the flight-path angle can be given by 


sin 7 = Z' 



( 8 ) 



6 


Thus^ the differential equation for the Z function during modulation 
contains the same terms as the eqixation developed in reference 1 for 
unmodulated entries plus the three additional terms in the curly brackets 
(cf. eqs. (4) and (6)). Equation (6) is applicable to both lifting and 
nonlifting entries. However^ as a result of the assimiptlons made in the 
derivation^ equation (6) is limited to shallow lifting entries for which 
I L/d tan 7 I « 1 (see appendix A) and is applicable to nonlifting entries 
for which the angles of descent are large as well as small. 

Equation (6) can be reduced to a differential equation in which u 
and Z are^ respectively, the independent and dependent variables if a 
drag polar and the type or purpose of the modulation are specified. In 
this manner it is possible to express L/D, A, A' , and A' ' in terms of 
Z, Z' , Z" , and u.^ Solution of the resulting equation for Z(u) permits 
the calculation of many quantities which describe the modulated portion 
of the trajectory; for example, resultant deceleration, flight-path angle, 
range, time, density, altitude, dynamic pressure, Reynolds number, con- 
vective heating rate, and total heat absorbed at the stagnation point. 
Detailed expressions for these quantities can be found in reference 1. 


Nonlifting Entries Modulated for Specified Maximum Deceleration 


As a specific application of the present analysis, modulated 
nonlifting entries will be studied. In this study m/CpA will be 
varied in such a manner that the maximum resultant deceleration does not 
exceed some specified value Gi. The modulated trajectory consists of 
three distinct phases which can be described with the aid of sketch (a) . 

The deceleration time histories of an unmodulated and a modulated entry 
are shoTO for a given entry angle (/p) and entry velocity (up) . During 
phase I the vehicle enters the atmosphere at a given entry angle and 
velocity (up) and maintains a constant m/CpA; the deceleration increases 
monotonically until the specified value of the deceleration Gx is attained 
at velocity Ux- At this velocity, phase II begins with the modulation of 
m/CpA in such a manner that the deceleration is held constant until the 
velocity is reduced to the value u^. This velocity is determined by 
specifying that termination of the modulation will allow the resultant 
deceleration to decrease monotonically as the vehicle coirpletes its entry. 
Phase III is this final phase of the entry. 


^or example, in an analysis to reduce the maximum resiuLtant 
deceleration Gmav , the equation for G given in reference 1 is used 
to relate G, L/d, Z, and u; that is, for shallow entries 

G = s/Fr uZ n/i + (l/d)^. 
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The present problem is reduced to the determination of the solution 
during the modulated portion of the trajectory (phase II) since the 
solutions for the unmodulated portions of the trajectory (phases I and III) 
are indicated in reference 1. There are three basic steps necessary to 
the solution for phase II. Firsts a relationship between Z, u, 7, 
and m/CpA must be obtained in order to determine the req^uired variation 
in m/CpA. Secondly^ the temlnation of the modulation period must be 
determined (i.e.^-the velocity Ua must be determined); and finally ^ the 
solution for phase II must be matched with the solutions for phases I 
and III. The first step is accomplished by means of the equation from 
reference 1 for the resultant deceleration; 

r = ^ 

cos^7 

During the modulation period G is constant and equal to Gi; hence 

^ uZ 

Gi r — 

cos-^7 


(9) 
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By combining eqioation (6) for L/D = 0 with equations (8) and (9) one 
can obtain a differential equation for the Z function in terms of pr_, 
U; and Gx^ but the mathematical treatment is tedious and the resulting 
equation is difficult to program for numerical solutions. This situation 
can be circumvented^ however^ if only shallow entries are considered; so 
that during phase II the flight-path angle is sufficiently small that 
cos^ ~ 1,^ ThuS; from equation (9) ^ the resultant deceleration during 
the modulation period for shallow entries is given by 


Gx 


(10) 


From eq-oatlon (lO) the Z fimction is defined during phase II as 



u 


( 11 ) 


■where 


B = 


Gi 


( 12 ) 


By substitution of equation (ll) and its derivatives into equation (6) 
(L/D = 0) , the differential equation which describes the variation in 
m/CpA required to maintain a constant resultant deceleration is obtained 
as 


A' ' 
A 


1 + A 

\Ay ■ -u ^ 


•w “ & 

~ A ■ COS^y = 0 


(13) 


and the related expression for the flight-path angle is 


sin 7 = 


B 

4 ^ 


"1 

/A'y 

2” 

u 

\A J 

"uA 




As a result of the multiplying factor cos^7 of the last term in 
equation (13) ; solutions for equation (I3) must be obtained by numerical 
methods . 

If d-oring phase II the flight-path angle is further assumed to be 
small enough that cos^7 « 1; a closed- form solution for A can be 


^Although the Inaccuracies due to assuming 003^7 « 1 were not 
investigated^ they are believed to be negligible for the entry angles 
considered in this report. 


-<] -fr* > 


4 =^ > 


2F 
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derived.. The resulting expression for the variation with velocity of the 
m/CpA required to maintain constant resultant deceleration is 



where 


C = 


hB^ 


> 


E = C(l - 2Gisin 7i - iii2) 


(16) 


The details of the derivation of equation (15) are presented in appendix B. 
The expression for the flight-path angle obtained with the closed-form 
solution is 


sin 7 = sin 7x + 


Gi 





(17) 


The accuracy of the closed-form solution relative to a numerical solution 
for equation (13) will he discussed later. 

The velocity ug at which the modiolation is terminated is determined 
in the following manner. For a given set of entry conditions and a 
specified deceleration limit, the constants B, C, and E are fixed, and 
it is a simple matter to determine from equation (15) that m/CpA 
increases to a maximum value and then decreases. Continuation of the 
modulation after m/CpA has reached the maximum would sustain a constant 
deceleration above that which would result if the modulation were termi- 
nated. Consequently, the value of the velocity U2 at which phase II is 
terminated is determined as that value for which A/a^ is a maximum. As 
noted in appendix B the value of U2 is obtained as that value of the 
velocity which satisfies the relation 





2C + 


B 


sin 7i 


= 0 


(18) 
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It now remains to matcli the solution for phase II wxth the solutions 
for the vinmodulated portions of the trajectory, phases I and III. This 
is done by matching the altitude and the flight-path angle at the veloci- 
ties for the beginning and end of modulation (ui and Ua, respectively) . 
The quantities Zi, or Ai, and 7x are known at velocity from the 

first phase of the trajectory (obtained as indicated in ref. l) ; hence, 
with the aid of equation (ll4-) the solutions for phases I and II are 
matched by satisfying the conditions 


^IIl 


= Ai 



4 ^ - . 2 
uxsin 7x + ~ 
B . ux 


> 


( 19 ) 


where the subscript Hi denotes initial values for phase II. The Z 
function and A are known at the end of phase II (eqs. (ll) and (l5)) and 
A'/a = 0 at velocity u^ (a/Ax is a maximum at u = ug) ; hence, with 
the aid of equations ( 5 ) and (l^) the solutions for phases II and III are 
matched by satisfying the conditions 


^Illi = ® 


Ug 


Z' III 


1 


( 20 ) 


which, together with the velocity u^, constitute the initial conditions 
from which the trajectory characteristics of phase III are calculated as 
indicated in reference 1. Within the basic approximations noted earlier, 
the closed-form solution is applicable to any vehicle making modulated, 
nonlifting, shallow entries into any exponential atmosphere from any 
approach velocity. 


Accuracy of the Closed-Form Solution 


The Inaccuracies due to assuming cos^y a 1 in obtaining the 
closed- form solution for phase II may be Judged from a con 5 >arison of 
some trajectory parameters for which the calculations for phase II were 
obtained from the closed-form solution and from the nmerical solution 
of equation ( 13 ) . The time history of the velocity u, the flight-path 
angle 7, and W/CpA^ are shown in figure 1 for a typical modulated. 


^Since the remainder of the report is concerned only with entries 
into the earth's atmosphere, the parameter W/CpA is subsequently used. 


^ > 
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lOg- limited^ nonlifting entry into the earth' s atmosphere at escape 
velocity (Y± = 1.4; 7 j^ = -8°) . Excellent agreement is exhibited between 
the solid curves representing the trajectory characteristics as computed 
from use of the numerical solution of equation ( 13 ) and the dashed curves 
representing the corresponding trajectory characteristics as confuted from 
use of the closed-form solution (eq. ( 15 )). 


RESULTS AND DISCUSSION 


The closed-form solution has been employed to study the effect of 
varying W/CdA on nonlifting shallow entries into the earth' s atmosphere 
A at escape velocity. In what follows, the effect of modulation on the 

4 time history of the resultant deceleration, flight-path angle, W/CpA, and 

4 velocity is discussed first. The effect of modulation on the corridor 

7 depth and heating during entry is then discussed. 


Trajectory Parameters 


Typical time histories of the resultant deceleration, flight-path 
angle, W/CpA, and velocity for nonlifting vehicles (l/D = O) are presented 
in figure 2 for unmodulated and modulated entries into the earth' s atmos- 
phere at escape velocity (Vi = 1.4) . The particular curves shown are for 
an entry angle - -6°. For the modulated entry, W/Cj)A was varied to 
maintain a peah deceleration of lOg duirlng phase II (see fig. 2(a)). It 
will be noted in figure 2(b) that the flight-path curvature during phase 
II of the modulated entry is considerably less than that for the unmodu- 
lated entry during the corresponding time interval. The reduced decelera- 
tion of the modulated entry results primarily from the reduced flight-path 
curvature. Since W/CpA is presented as a fraction of the value at the 
beginning of modiilatlon, the variation in W/CpA shown in figure 2(c) can 
be applied to vehicles of any weight, size, and shape. The relative areas 
under the velocity time-history curves (fig. 2(d)) are indicative of the 
greater range for the modulated entry. 


Corridor Depth 


In reference 2 the corridor depth in statute miles Ayp for 
successful entry into an exponential earth atmosphere was shown to be 


Ayp = 10 logic 


¥ 


P CpA 


un 


F- 


W 


P CpA 


ov 


( 21 ) 
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where the subscripts un and. ov , respectively^ 
and overshoot houndaries of the corridor and 
parameter given by 


Fp 


ui ^ 


refer to the "undershoot 
Fp is the conic perigee 


( 22 ) 


for 

_ ri _ 1 + .71 - Vi ^(2 - Vi^)cos^7i 

^ ~ ” -Tt- 2 2 

^ cos 

For the case of a nonlifting vehicle for which the initial 
the same for both boundaries 




10 logio 


^Pun 

^Pov 


W/CpA 


( 23 ) 

is 

(2h) 


and for a given entry altitude and vehicle (Zj^) > and a given entry velocity y 

(ui = ViCos 7 -) the corridor depth depends only upon the allowable entry 
angle (see eql. {22), { 23 ), and (24)), Thus the Increased corridor depths 
subsecjuently presented re suit from lowering the undershoot boundary by 
using modulation which permits deceleration-limited entries at steeper 
angles than is possible without modulation. 

Corridor depths and the corresponding entry angles attainable by 
modulating W/CpA during shallow nonlifting ^ entries into the earth's 
atmosphere at escape velocity are shown in figure 3 as a function of the 
ratio of w/CpA at the end of modulation to W/CpA at the beginning of 
modulation. Results are shown for deceleration limits of 15s* 

It can be seen that for a given deceleration limits the larger the avail- 
able change in W/CpA^ the broader is the corridor into which entry may 
be mad.e . For example, the depth of the lOg- limited corridor for^ 

an "unmodulated entry is Increased to 3 O miles for a modulated entry during 
which W/CpA at the end and beginning of modulation changes by a factor 
of 21. It is interesting to note from reference 2 that in order to use 
lift to attain a lOg-limited, 30-mlle corridor without modiilation, the 
entering vehicle must have an L/D of +O. 3 . It is also shovn in refer- _ 
ence 2 that without modulation the minlmimi peak deceleration for ^ nonlifting 
entry into the earth's atmosphere at escape velocity is 6 . 6 g, which occui s 
for 7 i = -4°. Consequently, no calculations for a 5g corridor have been 

made for entry angles less than -4°, 


Heating 


For the entries considered herein, peak heating occurs during the 
high-altitude portion of the trajectory where the densi"fcy and^ attendant 
Reynolds nnembers are low; consequently, the following discussion of 
heating d"uring entry will cover only those res"ults obtained when a laminar 
boundary layer is considered. 

Convective heating rate at th e stagnation point .- The convective 
heating rate per uniTlirea at the stagnation point during entry into the 


4^ 4T- > 
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h 

1 


earth's atmosphere ^ as given In reference is 

q, = 590 rn 

vhere 

for phases I and III 

5 = ' 

for phase II 


(25) 


( 26 ) 


Eq-uation (25) has Been used to calculate the stagnation-point heat-transfer 
rate for modulated and unmodulated entries which have the same ma:x;imum 
deceleration. Res-olts for the maxira'um heating rate are shown in figure 4 
for deceleration limits of 10 and 15g* In making these calculations it 
was assumed that the radius of curvature R of the vehicle surface at the 

stagnation point was identical and constant for hoth the modulated and 

unmodulated entries. It was also assumed that the initial W/CpA was the 
same for the modulated and the unmodulated entries. It should he noted 
that the data presented in figure 4 can he used to calculate the ratio of 
stagnation-point heating rates for the same modulated and unmodulated 
entries for which Rj^od / ^unmod^ provided R remains constant during 
entry. This is done hy multiplying the ratio of heating rates shown in 
figure 4 hy the square root of the ratio of the radius of curvature of 

the vehicle making the unmodulated entry to that of the vehicle making the 

modulated entry. If R varies during entry^ the stagnation-point heating 
rate will depend upon the variation of R with time or velocity (see 
eqs. (25) and (26)). 


The heating- rate penalty incurred hy the use of modulation (to 
increase corridor depth) can he seen in figiore 5 which shows the ratio of 
maximum heating rates for modulated and unmodulated entries as a function 
of corridor depth. For example^ in the case cited previously ^ in which 
the depth of a lOg-limited corridor is increased from 7 ho 30 miles ^ the 
maximum heating rate for the modulated entry which provides the 
corridor is 90 percent greater than that for the unmodulated entry which 
provides the 7-inile corridor. 


Total convective heat ahsorhed at the stagnation point .- The total 
convective heat ahsorhed per unit area at the stagnation point during a 
modulated entry into the earth^ s atmosphere can he written (see ref. l) 


where 


Qs = ±2^ f J 


= 19 ,900 
4^ 


W fdQ 


CdAE \du^ 


du 


(27) 


du 


r,3/2 


u 


„i/2 a 

Z cos 7 


for phases I and III 


B 


1/2 


(28) 


for phase II 



Equation (27) has "been used to calculate the total heat absorbed per unit 
area at the vehicle stagnation point for modulated and unmodulated entries 
which have the same deceleration limit. The results are shown in figure 6 
for deceleration limits of 10 and 15g- In the calculations the same 
assumptions were made with regard to vehicle shape and W/CpA that were 
made in calculating the heating- rate data of figure 4 (i.e., f^od = ^-unmod 
and (w/CdA)x = (w/Cj)A)xinmod) • vehicles which satisfy these assump- 

tions^ the penalty in total heat absorbed at the stagnation point can be 
seen in figure 6 . As in the case of the heat-transfer rate, the data of 
figure 6 may be similarly used to calc'ulate the total heat absorbed at 
the stagnation point for different vehicles when I^od ~ constant ^ Runmod^ 
otherwise R(u) must be available for use in equation ( 27 ) . 


CONCLUDING .REMARKS 


The differential equation developed for entries during which the 
aerodynamic coefficients and vehicle shape vary applies to vehicles of 
arbitrary weight, size, and shape and to arbitrary atmospheres. The 
closed-form solution for the motion, heating, and variation of W/CpA 
was obtained by making certain approximations which do not compromise the 
accuracy of the results. It was used to calculate trajectory parameters 
of nonlifting entries into the earth's atmosphere at escape velocity. As 
an example of the results obtained, the 7-mile lOg-llmited corridor for 
an unmodulated entry can be increased to 30 miles for a modulated entry 
with a W/CdA change of 21. For this same increase in corridor depth, 
however, the maximiim, laminar, convective heating rate and the total heat 
absorbed at the stagnation point are increased by factors of I .9 and 1 . 8 , 
respectively. These heating results were calculated for the case in which 
the radii of curvature were constant and equal during the modulated and 
unmodulated entries, and the initial values of W/CpA were identical. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif, Aug. 12, I 96 O 
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APPENDIX A 


DIFFERENTIAL EQUATION OF THE Z FUNCTION FOR MODULATED ENTRIES 


The development of the differential equation of the Z function^ 
for modulated, tvo -dimensional (no lateral forces) entries into spherically 
symmetric, exponential, planetary atmospheres closely parallels that pre- 
sented in reference 1 for unmodulated entries. Deviations from the 
analysis of reference 1, it "will he seen, occur only vhen velocity deriv- 
atives are obtained, since, in the present analysis, the aerodynamic 
forces and vehicle shape (and weight) are considered to depend on velocity. 

A. 

4 By use of a polar coordinate system with velocity components, 

4 aerodynamic forces, and flight-path angle defined as in sketch (h) , the 

7 



pair of motion equations in the radial and tangential directions, 
respectively, can he shown to he 



d^y _ dv 


L 

COS 

D . 

y + — sin 7 

(Al) 


" dt^ “ " dt 

® r 

m 


' m 



du ^ uv _ 

D / 

- cos 7 ( 

1 


tan 7 ) 

(A2) 


dt r 

m ^ 


D 

J 



By utilizing the same approximations used in reference 1 so that uv/r can 
he neglected compared to du/dt (i.e., | dr/rl « Idu/ul)^ hy Introducing 
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the drag coefficient D = (l/ 2 )Cj)pV^A and the exponential atmosphere 
(p = Po®”^^) ? t’y noting that V = u/cos y, and by introducing the inde- 
pendent variable 


u = 


u 




(A3) 


so that 


du 


1 du 


dt 4 ^ dt 

one can write equations (A 1 ) and (A 2 ) 


_ i ^ ^ = 1 _ 


g dt" 


g dt 


+ 


rp^e 


m 


CdA 


cos^7 


sin 7 - 


- cos 7 


(aU) 


(A5) 


du j P Pn^ 


dt= _ /j5_\ 




1 + ^ tan 7) 


(A 6 ) 


In order to reduce the pair of motion equations (A5) and (a 6 ) to a 
single equation^ the dependent variable Z defined by 


1 _ fr e"^^u 

2 Po J 13 ^ 

CdA 


(A?) 


is introduced. With the aid of equation (AY) equations (A5) and (A 6 ) are 
more conveniently written 


1 ^ 
g dt 


1 - u^ + 


4 Jv uZ 
cos^7 


sin 7 


— cos 
D 



(A8) 


dt 


cos 7 \ D y 


(A9) 




-<l -p- -r > 


-<1 -p -r > 


3F 
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A second equation for dv/dt is developed in terms of Z, j , and 
m/CpA by the use of equations (Ay), i^ 9 ) > and the following expression 
for the radial component of velocity (see sketch (b)) 


iz = 

dt 


V = gr u tan 7 


(AlO) 


The resulting equation is then equated to equation (A8) to obtain the 
desired single equation in terms of the foregoing variables. 

Before the second equation for dv/dt is developed it is convenient 
to introduce the notation 


A - 


m 

CqA 


so that 


A* 



du d1I 


Thus^ equation (AY) is rewritten 


Z = 


1 

2 


Po 


/7 e-PZu 

n/ P A 


(All) 


(A 12 ) 


(AI 3 ) 


Proceeding now with the development, differentiation of equation (AlO) 
yields 


]z - [^ d(u tan 7) _ fi. ^ f- sin du 
g dt J g dt ~ J g d.u\^ cos yj dt 


(Al4) 


Upon substituting equation (A 9 ) in equation (A 14 ) and noting that 

^ ^ ^ , we obtain 

du dy du 


1 ^ 
g It 


uZ 

cos^7 


(n/Pt 


V 


sin 7 + 


>/pr u d7\ 
cos 7 diy 



(AI5) 


To provide the desired foirni of the second expression for dv/dt (eq. (AI5) ) , 
the flight-path curvature d7/du is expressed in terms of u, Z, 7, and 
A in the following manner. By differentiating Z/u from equation (AI3) 

(Z' = dZ/dtl) and keeping in mind the approximation 1 dr/r| « | du/^l 
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Z' Z _ 1 - fr £fy f. ^ 

u " " 2 J P A V du A, 


Z ^0 ^ M + ^' 

-"- ^ ^ dt da A 


u 


(Al 6 ) 


Use of equations (A 9 ) and (AlO) in equation (Al 6 ) yields an expression 
for the flight -path angle 


sin 7 = ( Z' - I + Z (l +|tan 7 


(AI7) 


An expression for the flight-path curvature d7/duV obtained by differen- 
tiating equation (AI7) , is .. 


A 

4 

4 

7 ’ 


1 + — tein 7 j u ^ 
D V du 


U 


n/^ cos 7 




7. Z ^ COS.27 

^ "u”^ A J ypy cos 7 


- 'i 


tan 7 


Xj \ 2 

1 -h — tan 7 ) cos 7 


(AI8) 


as Sliming that L/d tan 7 can he neglected compared to unity ^ 

(1 L/D tan 7i « l) and restricting the analysis for lifting vehicles to 
shallov entries so that during the modulation period cos^7 - 1; ve can 
reduce equation (A18) to 

u ^ = _ — H ^fz' - I + z ^ 

du „ypr cos 7 \ u Ay 



Finally, upon substituting equation f AI9) in equation (AI5) (l/D tan 7 « l) 
and equating the result to equation (a 8) the original pair of motion 
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equations (Al) and (A2) are transformed into a single differential equation 
for modulated entries into exponential planetary atmospheres 


tiZ' ' 



A' ' 

A 



+ uZ' 



A 

4 

7 


1 - u2 


uZ 


COS^y + — cos^7 = 0 


(A20) 


In this equation; cos 7 = ^1 - sin ^7 can he expressed in terms of Z, 
Z' ; A; and A' through equation (AI 7 ) (l/D tan 7 « l) . 
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APPENDIX B 

CLOSED-FORM SOLUTIONS FOR m/CpA, ua, AND 7 FOR 
MODULATED NONLIFTING ENTRIES 

The variation in m/CoA with velocity required to maintain a 
constant resultant deceleration is defined hy equation (13) • Upon suh- 
stitution of the dependent variable 


_ ^ ^ 1 ^ ^ 

P “ A A du “ m/Ci)A 


d(m/Cr)A) 


du 


and assuming that cos^y 1, equation ( 13 ) reduces to 


(Bl) 


p' - - P = 


1 - u^ 


B" 


u2 


(B2) 


The corresponding expression for the flight-path angle is (see eqs. (Bl) 
and ( 1^) ) 


'JJr sin 7 = |(p - f 


(B3) 


Integration of equation (B2) yields (l/u is an integrating factor) 


i p = — in u - ~ + ^ + constant 


u 


B" 


2B^ u^ 


(b4) 


The integration constant can he evaluated in terms of the velocity 

the flight-path angle at the beginning of modulation (i-e., from eq. 1, j,) , 

at u = ui, Pi = (Wb) ^ sin 71 + ( 2 /Ui)) - In this maimer 


A' _ _ n , ^ f i 

-7 P " ~ ^ n ( - _ j + ^^2 V "1 1 




Ul 


1 - ( 

ni 


+ I + u sin 7i (B 5 ) 
n B 


Proceeding ty iniegra'tion of ec][ua'bion (Bl) 


^ ^ = exp y jp du -H constanij 


(b6) 


-j -r > 
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A 

4 

4 

7 


Equation (B 5 ) is readily integrable and the integration constant can be 
evaluated in terms of the velocity and m/Ci)A at the beginning of 
modulation (ui and Ai, respectively). Thus, the closed-form solution 
for the variation of m/CpA during a nonlifting modulated entry is 


^ m/CpA 

(hi/CdA) 1 


exp 



+ I CUi2 


1 



(B7) 


where 


C 


4B^ 


E = C(l - 2 B sfpr sin 71 - Ui^ 


(b8) 


and B is given by equation (l 2 ) . 

The velocity u.2 at which the modulation should be terminated has 
been shown to be that value at which m/Cj)k reaches a maximum, that is, 
the value of u for which the derivative of equation (b 6) (eq. (B5)) is 
zero, which for convenience is written 


B^ 




(b9) 


Trial and error solutions of equation (B 9 ) yield the velocity U2- The 
expression for the flight-path angle during the modulation period is 
obtained from substitution of equation (B5) into equation (B3) 


sin 7 = sin 71 





1 



(BIO) 
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(c) Modulation parameter 





Figure 2,- Concluded 
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